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Abstract. In this paper we study, both analytically and numerically, questions involving the distribution of eigen- 
values of Maass forms on the moonshine groups To(N) + , where N > 1 is a square- free integer. After we prove 
that Fq(N) + has one cusp, we compute the constant term of the associated non-holomorphic Eisenstein series. We 
then derive an "average" Weyl's law for the distribution of eigenvalues of Maass forms, from which we prove the 
"classical" Weyl's law as a special case. The groups corresponding to N = 5 and N = 6 have the same signature; 
however, our analysis shows that, asymptotically, there are infinitely more cusp forms for To(5) + than for To(6) + . 
We view this result as being consistent with the Phillips-Sarnak philosophy since we have shown, unconditionally, 
the existence of two groups which have different Weyl's laws. In addition, we employ Hejhal's algorithm, together 
with recently developed refinements from | 31| . and numerically determine the first 3557 of To(5) + and the first 12474 
eigenvalues of ro(6) + . With this information, we empirically verify some conjectured distributional properties of 
the eigenvalues. 



Let {pi}, with i = 1, . . . , r, be a set of distinct primes, so then N = p\ ■ ■ ■ p r is a square-free, non-negative integer. 
The subset of SL(2,M), defined by 



is an arithmetic subgroup of SL(2, M). The groups To(N) + were first considered by Helling [TS] where it was proved 
that if a subgroup G C SL(2,M) is commensurable with SL(2,Z), then there exists a square-free, non-negative 
integer N such that G is a subgroup of Tq(N) + . We also refer to page 27 of [27] where the groups Tq(N) + are 
cited as examples of groups which are commensurable with SL(2,Z) but non necessarily conjugate to a subgroup 



Following the discussion in [TU] [HJ HH E3 E2 , we employ the term "moonshine group" when discussing To(N) + . 
The genus zero moonshine subgroups of SL(2,R) arise in the "monstrous moonshine" conjectures of Conway and 
Norton, which were later proved in the celebrated work of Borcherds. Gannon's book [15] provides an excellent 
discussion of the mathematics and mathematical history of monstrous moonshine. In particular, we refer to Con- 
jecture 7.1.1 where the Conway-Norton conjecture is stated, which in its original form referred to certain genus zero 
subgroups of Moonshine-type. After the work of Borcherds, the authors in [TU] described solely in group-theoretic 
terms the 171 genus zero subgroups that appear in mathematics of "monstrous moonshine" . Amongst this list are 
those groups of the form ro(-/V) + which have genus zero. 

Our interest in the groups To(iV) + stems from the work in [23]. In that article, the groups ro(5) + and To(6) + 
were examples of arithmetically defined topologically equivalent groups which have distinct spectral properties. 
More specifically, in |23j the authors defined an invariant associated to any non-compact, finite volume hyperbolic 
Ricmann surface, where the invariant is equal to the larger of two quantities: one coming from the length spectrum 
and another associated to the determinant of the scattering matrix. The groups ro(5) + and To(6) + have the same 
signature and are arithmetically defined, yet have different values of the invariant defined in |23j . As a result, the 
main theorem of 23| showed that, in somewhat vague terms, the derivative of the Selberg zeta function of one 
surface has more zeros than the derivative of the Selberg zeta function of the other. Since the spectrum of a surface 
is measured by the zeros of the Selberg zeta function, the main result of [5H| can be interpreted as saying that 
surfaces r (5) + and r (6) + are quite different from the point of view of the asymptotics of spectral analysis. 

In somewhat vague terms, the purpose of the present article is to investigate the spectral properties of the 
Ricmann surfaces associated to the groups r (A^) + for square- free N in order to make precise the observations 
made in 23 . In doing so, we employ the ideas from |31j which build on Hejhal's algorithm for numerically 
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estimating eigenvalues of the Laplacian on finite volume, hyperbolic Riemann surfaces. With this said, we now can 
describe the main results. 

Let To(N) + — To(N) + / ± /, where I is the identity matrix and let Xn :— Yq(N) + \B. be the corresponding two 
dimensional surface. Since T (N) C T (N) + , where T (iV) denotes the classical congruence subgroup of SL(2,Z), 
the surface Xjy has finite volume. As stated, we will show that for any square-free TV, the surface has exactly 
one cusp; hence the signature of r (7V) + is (g; m\, . . . ,ttij; 1) where g denotes the genus of the group and I is the 
number of inequivalent elliptic elements of To(iV) + with mj, i — 1, . . . , I denoting the order of the corresponding 
elliptic element. 

Maass forms on Tq(N) + are real analytic, square integrable, eigenfunctions of the Laplacian on the surface Xn. 
Maass forms which vanish in the cusp are called Maass cusp forms. The hyperbolic Laplacian —A on Xn has a 
discrete and continuous spectrum; see |21) or |17| . The discrete spectrum is denoted by the set {A„} n >o, counted 
with multiplicities; here, we have that = Ao < Ai < . . . < A„ N _i < 1/4 < A, iN < . . . and A„ — > 00 as n — > 00. Let 
m i/4.N > denote the multiplicity of A = 1/4 as (eventual) eigenvalue of —A. Maass cusp forms span the positive 
discrete part of the spectrum. 

Let {r n } denote the set of all positive real numbers satisfying the equation 1/4 + r\ = A„. For T > 0, the 
function ~Nn(T) := Njv[0 < r n < T] counts the number of r n such that < r n < T, or, equivalently, the number of 
eigenvalues of Maass cusp forms which lie in the interval (1/4, T 2 + 1/4]. 

For any T > and square-free N, which we write as N = p\ • • -p r , define ajv(j, T) :— T log pj — [ T — Pj j n where 
[x\ denotes the greatest integer less than or equal to x. 

The main analytical result of the paper is the following theorem: 

Theorem 1 (Average Weyl's law for r n (iV) + ). Let (g; mi, . . . , m;; 1) be the signature of the group r (-/V) + and let 
nj\r > 1 denote the number of small eigenvalues of the Laplacian — A on Ajy. Then 



N N (T)=M N (T)+S N (T) 



where 



M N (T) = M^T^ - 27 ^ gT + T (2 + l 0g(7r /2iV)) + V JL £ MM 

47r 7r 7T f-^ 4m t ^— { sin (nj/mi) 487T 

3 ri]\[ 1 \ - . 1 x ^ 

_ 4 ~ + 2*2^ aN{h T) - ^2^ arctan 

3=1 i=i 




1/4, N 



G N (T), 



with 



1 / Vol(A jv )(2^ + l) ' m t m ^ 1 5051 \ 1 

U 1 N[ ,] ~ 2tt I 27r 2 cxp(27r) ^ 2 e7 r ^ si^Trj'/m,) 900 J T' 

/or all T > 1 and 

T 

j S N (t)dt = ( T 2 J j as 00. 


The word "average" in the title of our main theorem relates to the form of the error term in the Weyl's law. 
An average Weyl's law is of importance when it comes to the numerical computation of Maass forms; see [31] and 
references therein. In particular, when computing Maass cusp forms numerically, there is always the risk that some 
solutions get overlooked. By comparing a numerically found list of eigenvalues of Maass cusp forms with average 
Weyl's law, one can easily determine the number of solutions which have been overlooked. We refer to |31| for a 
detailed discussion of this point. 

An immediate consequence of Theorem [l] and its proof is the following corollary. 

Corollary 2 (Classical Weyl's law for T (N)+). 

K„m . - + T -(2 + ■„ gW 2iV» + O . « T -, 00. 

47T 7T 7T \lOgl ) 

Generally speaking, the philosophy behind the Phillips-Sarnak conjecture [25j [26] suggests that the spectral 
analysis of the Laplacian acting on smooth functions on a finite volume, hyperbolic Riemann surface M should 
depend on the arithmetic nature of the underlying Fuchsian group T. The first terms in the asymptotic expansion 
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in Corollary [2] depend solely on the volume of Xjv, and then one sees that the coefficient of T depends on N. For 
example, the groups corresponding to N = 5 and N = 6 have the same signature, hence X§ and Xq have the same 
volume yet, by Corollary [2j X§ has infinitely more eigenvalues than X§ in the sense that 

lim I (X 5 (T) - N 6 (T)) = log(6/5) > 0. 

Later in this article, we provide a list of further examples of topologically equivalent surfaces associated to moonshine 
groups which have different Weyl's laws. We view these results as being consistent with and in support of the 
Phillips-Sarnak philosophy. 

Having established that the classical Weyl's law associated to ro(5) + and To(6) + differ, we find it interesting 
to investigate other conjectures concerning the distribution of eigenvalues. Using the methodology from |31j . and 
references therein, we have numerically computed sets of Maass cusp forms associated to ro(5) + and To(6) + . On 
To(5) + our numerical results cover the range < A < 125 2 + 1/4 which includes 3557 Maass cusp forms, and on 
r (6) + we cover the range < A < 230 2 + 1/4 which includes 12474 Maass cusp forms. The distribution of the 
numerically found eigenvalues is in agreement with the following conjecture. 

Conjecture 3 (Arithmetic Quantum Chaos [31 |S]). On surfaces of constant negative curvature that are generated 
by arithmetic fundamental groups, the distribution of the discrete eigenvalues of the hyperbolic Laplacian approaches 
a Poisson distribution as A — > 00. 

A particular feature of a Poisson distribution is the "absence of memory" , which, in our case, asserts that an 
eigenvalue cannot be predicted from knowledge of all the previous eigenvalues. The computation of eigenvalues 
allows us to verify that, numerically, eigenvalues of the Laplacian on A5 and X$ are uncorrelated. 

This paper is organized as follows. In section [2] we provide preliminary material for both the theoretical and 
numerical aspects of our work. Theoretically, we prove that the Ricmann surfaces associated to the moonshine 
groups To(A^) + for square-free N have one cusp, and we compute the first Fourier coefficient of the corresponding 
non-holomorphic Eisenstein series. In order to make this article as self-contained as possible, we include a discussion 
of Hejhal's algorithm for numerically estimating eigenvalues together with Turing's method which is used to verify 
that no eigenvalue has been missed. In section [3] we prove Theorem [l] and as corollaries state the result in the 
cases of SL(2,Z), To(5) + and To(6) + . In sectionHwe state the conclusions from our numerical investigations, and 
in section [5] we present various concluding remarks. 



2. Preliminaries 

2.1. Moonshine groups T (N) + . In this subsection we will derive some important properties of moonshine groups 
r (-/V) + , for a square-free integer N. We will prove they have exactly one cusp. We then compute the constant 
Fourier coefficient of the associated non-holomorphic Eisenstein series. Equivalently, we compute the scattering 
determinant associated to the cusp. We refer to |17j and |21j for relevant background information. 

Lemma 4. For every square-free integer N > 1, the surface Xn has exactly one cusp, which can be taken to be at 
ioo. 

Proof. The cusps of Xjy are uniquely determined by parabolic elements of the group r (Af) + . In [T2] it is proved 
that all parabolic elements of T${N) + have integral entries. Therefore, the parabolic elements of T${N) + are also 
parabolic elements of the congruence group T^(N). From pages 44-47 of J2T], we easily deduce that the only possible 
cusps of ro(A^)+\H belong to the set {0, zoo} U {1/v : v \ N}. The point z = is mapped to zoo by involution 



For an arbitrary v \ N and w — N/v one has (w, v) = 1 since N is square-free. By Euclid's algorithm, there 
exists integers a and b such that —aw — bv = l. Therefore, points z = 1/v are mapped to zoo by transformation 



1 / aw b 

~ \ N -w 



€ T (N)- 



This shows that all possible cusps of ro(iV)+\H are To(iV) + -equivalent with zoo. Therefore Fo(iV)+\H has exactly 
one cusp which can be taken to be zoo, as claimed. □ 



Let C(s) denote the (classical) Riemann zeta function and let £(s) be the completed zeta function, defined by 
Z(s):=is(s-l)ir-*/ 2 r( S /2)(( S ). 
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Lemma 5. For a square-free, positive integer N = p%- ■ -p r , the scattering determinant associated to the cusp of 
Xn at ioo is given by the following expression 

, % s £(2s — 1) „ , , 

(2) = 7 c (0 s ■ Dn ^> 

s — 1 s(2sj 

where 

Proof. By Theorem 3.4 from [3T] we write <pn{js) — ^JttT(s — \/2)T~ 1 (s)Hn{s) 1 where fljv(s) denotes the Dirichlet 
series portion of the scattering determinant. Let Cat denote the set of left-lower entries of matrices from Tq(N) + . 
Following pages 45-49 from |21j . one sees that 

Hn(s) = c ~ 2s An(c) 

c£C N 

is well defined for Re(s) > 1, where Am{c) is equal to the number of distinct values of d modulo c such that c and 
d are elements of the bottom row of the matrix from Tq(N) + . 

From the definition of T Q (N) + , we easily deduce that Cjv = {(N/^/v) ■ n : v \ N,n € N}. 

For a fixed c = (N/y/v) ■ n, with v | N and n € N arbitrary, we can take e = u in the definition of r (iV) + to 
deduce that matrices from Tq(N) + with left lower entry c are given by 

y/va b/y/v 
^y/vn y/vd 

for some integers a, b and d such that vad — (N/v)bn = 1. Therefore, the number An((N/-^/v) ■ n) is equal to the 
number of distinct solutions d modulo {N/v)n of the equation vad — (N/v)bn = 1. Since N is square-free, this 
equation has a solution if and only if (v,n) = 1 and (d, (N/v)n) = 1. In this case, the number of distinct solutions 
d modulo (N/v)n is equal to <p((N/v)n). Here, <p denotes the Euler totient function and (p,q) denotes the greatest 
common divisor of integers p and q. 

Therefore, A N ((N/y/v) ■ n) = if (v,n) ^ 1 and A N ({N/y/v) ■ n) = <p((N/v)n) if (v,n) = 1. Now, we may 
conclude that 



u|JV (n,i;)=l (irV" Tl ) 



2 s ' 



The inner sum on the right-hand side of the above equation may be expressed using computations from [17] . 
specifically Lemmata 4.5 and 4.6 on page 535, showing that for positive integers A\, A 2) B\ and B2 one has 



¥ >{c -(B 1 ,B 2 )(A 1 ,A 2 )) 
c 2 n s (A ll A 2 ) s (B 1 ,B 2 ) 2s 

c >0: (c ,{B 1 ,A 2 )(A 1 .B 2 )) = l V U 2; V 1, M 



(3) E 



c(2 - 1) n (^) n ^'- p '" 



C(2s) 11 VP 25-1 J 11 V P 2s_1 

in standard notation, p denotes a prime number, and an empty product is defined to be equal to 1. 

Using formula ([3| with A\ = v, A 2 = 1; B\ = iV/u, B 2 = N and the principle of mathematical induction with 
respect to the number r of distinct prime factors of N = p\ • • -p r , we deduce that 

C(2«-l) V [ j-r P-I jt P'-P 1 -" j = C(2«-1) 1 fr Pl+Pj 
^ W C(2s) ^ 11 p2«-lll a( 2s ) llpj + l' 

s ^ 7 "1^ \p|(f ) pI« / J 

Therefore, the scattering matrix y>jv(s)j f° r Re(s) > 1 is given by 

_r(s-l/2) C(2s- 1) 

r(s) c(2s) 

The statement of the lemma follows from the definition of the completed zeta function, which completes the proof 
of the Lemma. □ 

Remark 6. The determinant of the scattering matrix for congruence subgroups has been computed by Hejhal [T7] 
and Huxley [H]. 
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Figure 1. Dirichlet fundamental domains of the moonshine groups ro(5) + (left), and To(6) + (right). 



2.2. Moonshine groups T (5) + and r (6) + . The moonshine group r (5) + is generated by 

' \ 1\ I /'•-. -1\ 1 ( -:-! 

and the moonshine group To(6) + is generated by 



1 l\ 1/6 -1\ 1 /3 -2 



see [T3]. Fundamental domains of X$ — Fo(5)+\H and Xq — ro(6) + \H are displayed in figure]!] For both, X$ and 
Ag, the sides are identified according to the pairings 

gx : si s 6 , 32 : s 2 H> s 5 , g 3 : s 3 ^ s 4 . 

Both A5 and Ag have a cusp at «i = ioo, and each surface has three inequivalent elliptic fixed points which are all 
of order 2. The elliptic fixed points are 

9\ 1 Q2 '■ V2 1— > V2 which is T-equivalent with vq = g\V2, 
92 X 9z '■ v 3 ^ v 3 which is T-equivalent with v$ — .93 ^3, 
and 33 : V4 W4. 

By the Gauss-Bonnet theorem, the volumes of the surfaces are Vol(X5) = ir and Vo\(Xq) = it. 

2.3. Strombergsson's pullback algorithm. In 2000, Strombergsson [29] presented an algorithm for computing 
the pullback of any point z G H into the Dirichlet fundamental domain of a given cofinite Fuchsian group T with 
prescribed generators. Strombergsson's algorithm uses only the action of generators of the group T applied to the 
point z, and the algorithm is shown to converge after a finite number of iterations. Computation of the pullback of 
a point z S H to the Dirichlet fundamental domain of F is an ingredient in Hejhal's algorithm for computing Maass 
forms, recalled below. Therefore, Strombergsson's algorithm is an important part of our numerical computations 
of eigenvalues of Maass forms on X$ and Xq. 

For the sake of completeness, we will recall the Strombergsson algorithm in its full generality. Assume that T 
is a cofinite Fuchsian group with generators g\,...,g n and set of elliptic fixed points £. Let d(z, w) denote the 
hyperbolic distance between two points z and w in H. The associated Dirichlet fundamental domain is the set 

J = d{p, z) < d{ P , jz) v 7 e r}, 

where p € HI — £ is arbitrary. The given generators of T identify the sides of 3\ Strombergsson's algorithm for 
computing the pullback of any point z € H into 7 is the following. 

Algorithm 7 (Pullback algorithm [HI]). Choose any z E EL 

(1) Compute the 2n points g±z, g^ 1 z, g2Z, g^ 1 z, . . . , g„ lz - Let z' be the one of these points which has the smallest 
hyperbolic distance to p. 

(2) If d(p, z') < d(p, z), then replace z by z' , and repeat with s£ep[7J 

(3) // d(p, z') > d(p, z), then we know that z lies in J ', hence z is the desired point, i.e. the pullback of the point 
initially selected. 

Strombergsson proved that his algorithm always finds the pullback within a finite number of operations [29] . 
We use z* — x* + iy* to denote the pullback of z = x + iy. 
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2.4. Maass forms on To(iV) + . Let us recall the definition of Maass forms [23] and Maass cusp forms. 

Definition 8. / : H — > R is a Maass form on ro(-ZV) + associated to the eigenvalue A if and only if 

i) /eC»(i), 

ii) / G L 2 (X N ), 

iii) -Af(z) = Xf(z), 

iv) fiiz) = f(z) V 7 G T (N)+. 

Definition 9. / : H — > R is a Maass cusp form on T (N) + if and only if 

i) f is a Maass form on To(N) + , 

ii) lim z ^ loo f(z) = 0. 

For z — x + iy G H, the Fourier expansion of a Maass cusp form associated to the eigenvalue A = r 2 + 1/4 is 
given by 

(4) f(x + iy)= ]T a^K^n^e 2 ^, 

nGZ-{0} 

where K stands for the if-Bessel function. Since a Maass form is real analytic, we have Rea_„ = Rea„ and 
Ima_„ = — Ima„. 

As first proved in [23] . the spectral coefficients a n grow at most polynomially in n. The if-Bessel function decays 
exponentially for large arguments, meaning 

K„{y) ~ \j^y e ~ V for y^°°- 

As a result, one can obtain a very good approximation of the expansion Q by using finitely many terms, where 
the number of terms considered depends on the desired accuracy of the approximation. 

Let 3 N ~ r o (A0+\H be the fundamental domain of T (N)+. Let z* = x* + iy* be the r (iV) + -pullback of the 
point z = x + iy into the fundamental domain, meaning there exists some 7 G r (-/V) + such that z* = jz and 
z* G Jn. By the definition of automorphy, we have that f(z) = f(z*). 

Since the congruence group Tq(N) is a subgroup of ro(-/V) + , we immediately deduce the following lemma. 

Lemma 10. If f is a Maass form on Tq(N) + , then f is a Maass form on Tq(N). 

2.5. Hecke operators. Let us recall the definition of Hecke operators. There are many references for this material, 
one of which being [27] . 

Definition 11. Let f : H — > R, and n a positive integer. The Hecke operator T n is defined by 

Jn z — ' L — ' 

ad—n 6—0 
d>0 

Theorem 12 ([Tll27j). Consider the congruence group Tq(N). For all n such that (n,N) = 1, the Hecke operators 
T n are endomorphisms of the space of Maass cusp forms on Tq(N). For all m and n with (m, N) — (n, N) = 1 and 
all Maass cusp forms f(z) on Tq(N), the Hecke operators have the following properties: 

T m T n — Tm^ , 

d\ (m,n) 

T n o A = A o T n , 

Tnf(z) = t n f(z), 

where the eigenvalues t n of the Hecke operators T n are related to the expansion coefficients a n of the Maass cusp 
form f{z) by the identity 

a n a\t n . 

For a proof of the theorem, see [1] , [27] , or [28] . 

Theorem 12 immediately implies that the Fourier coefficients of Maass cusp forms on T (N) are multiplicative, 



d - 1 , u 

,az + b. 



a 1 



d\ (m,n) 
d>0 



for all m and n with (m, N) — (n, N) = 1. By Lemma 10 this holds also for Maass cusp forms on T (N)~ 
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2.6. Hejhal's algorithm. We make use of Hejhal's algorithm [18, 30j which itself employs the Fourier expansion 
Q of Maass cusp forms. 

Hejhal's algorithm is a finite system of linear equations whose non-trivial solutions are related to Maass cusp 
forms. Hejhal's algorithm is heuristic. By construction, a Maass cusp form always will solve the linear equations of 
the algorithm to any desired level of accuracy, but the converse is not true. Not each solution of the finite system 
of linear equations is a Maass cusp form. Only in the case when a solution is independent of the parameters will 
the solution approximate a Maass cusp form. The crucial parameter in question is the choice of the value of y in 
Q. The computation of Maass cusp forms therefore proceeds in two steps: Heuristic use of Hejhal's algorithm, 
followed by a verification of the numerical results. 

Theoretically, Maass cusp forms can be rigorously certified as was shown in [9] in the example of the modular 
group. Using the quasi-mode construction, Booker, Strombergsson, and Venkatesh have certified the first 10 eigen- 
values of SL(2,Z). The certification techniques can be adopted to other settings, such as to the moonshine groups. 
Practically, however, we have to bear in mind that rigorously certifying eigenvalues requires immense computer 
resources and it is infeasible to certify thousands of Maass cusp forms. For this reason, we just verify the numerical 
results with a different, not fully rigorous method. 

The verification is based on the following: 

(1) Fix y. 

(2) Find non-trivial solutions of Hejhal's system of linear equations. 

(3) Take a finite number of different values of y, and check whether the non-trivial solutions seem to be 
independent of y. 

(4) Take only the solutions which are seemingly independent of y and make a list of conjectured Maass cusp 
forms. 

In the end, there will be strong evidence, but not a proof, that the list of conjectured Maass cusp forms is indeed 
a list of true Maass cusp forms. It is the experience of those who implement the algorithm that more than half of 
the non-trivial solutions of Hejhal's system of equations for a fixed value of y are not Maass cusp forms. Taking a 
second choice for y immediately rules out almost all solutions which are not a Maass cusp form. 

There remains the possibility that a solution could solve Hejhal's linear system of equations for two independent 
values of y whilest not being a Maass cusp form. We have further checked whether this has happened by employing 
several independent values of y. Empirically, it turned out that as soon as some function solves Hejhal's system of 
equations for two independent values of y, it does so for any finite number of independent values of y also. And we 
conjecture that it does so for any other value of y. 

Further evidence comes from a second verification based on the Hecke operators. According to the Hecke 
operators, the expansion coefficients of Maass forms are multiplicative. When solving Hejhal's system of linear 
equations, there is no reason that the coefficients of a solution are multiplicative, but only those solutions whose 
coefficients are multiplicative can be Maass cusp forms. 

Numerically, for each individual solution of Hejhal's system of linear equations we have investigated and found 
that a solution is seemingly independent of y if and only if the expansion coefficients of the solution are multiplicative. 
This means both verifications agree in their answer. 

Let us now recall Hejhal's algorithm. 

Since Tq(N) + is cofinite and has only one cusp at zoo, we can bound y from below. Allowing for a small numerical 
error of at most [[e]], where [[e]] stands for \numerical error\ < e, due to the exponential decay of the K— Besscl 
function in y, we can truncate the absolutely convergent Fourier expansion Q such that 

(5) f(x + iy)= J2 a n y l / 2 K ir {2ir\n\y)e 2 * inx + M\. 

0jt\n\<M(e,r,y) 

Solving for the spectral coefficients results in the equation 

, 2Q 

(6) a m y^K tr (2n\m\y) = — £ f{j- + ? y) e - 2 ™"* + [[e]}, 

3=1 

with 2Q > M + to. 

By automorphy, any Maass cusp form can be approximated by 

(7) f(x + iy) = f(x*+iy*)= £ a n y* 1/2 K lr {2^\n\y*)e 2mnx * + [[e]], 

where y* is always larger than or equal to the height of the lowest point of the fundamental domain 5F, allowing us 
to replace M(e,r,y*) by M = M(s, r, min^ggTmiu). 
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Making use of the implicit automorphy by replacing f(x + iy) in ^ with the right-hand side of ^ yields 

(8) 

a m y X/2 K ir {Mm\y) = — ^ E a n y* 1/2 K tr (2n\n\y*)e 2 ^ nx *i- m ^ + [[2e]], where Xj + iy, = j- + iy, 

for ^ |m| < M(e, r, y), which is the central identity of the algorithm. 

We are looking for non-trivial solutions numerically such that ^ vanishes simultaneously for all 7^ \m\ < Mq 
and < y < min^gjlmw. Each non-trivial solution gives a Maass cusp form whose eigenvalue reads A = r 2 + 1/4. 

We first solve Q for all 7^ \m\ < Mq numerically, but use a single value of y only. Then, we verify with a 
finite number of values of y, whether we have found a non-trivial solution such that ([8]) vanishes simultaneously 
for all 7^ \m\ < Mo for each value of y. If the solution turns out to be seemingly independent of y, we finally 
check whether the expansion coefficients a n are multiplicative. If also the expansion coefficients turn out to be 
multiplicative, we have verified that the numerically found solution of ^ is a Maass cusp form. 

Let us now specify good parameter values for solving ^ numerically. 

Algorithm 13 (Parameter values). Let A = t 2 + 1/4 be close to an eigenvalue. Let the precision be given by e > 0. 
Then for A near A we choose the values of the parameters as follows: 

(1) Solve eK it (max{t, 1}) = K it (2irM Q min^ggrlmw) in Mq with 2TrMomin w ^hxiw > max{t,l}. 

, 9 maxjt, 1} 

2 Let y = 

w y 10 2-kMq 

(3) Solve eK it (max{t, 1}) = K it (2nMy) in M with 2-nMy > max{i, 1}, i.e. M = mm ™£3 rlmw Mq _ 

y 

(4) Let Q be the smallest integer which is larger than M . 

(5) Check whether ^ is well conditioned for the given y and all =/= \m\ < Mq. Lf not, reduce y slightly and 
repeat with^ 

For verifying that ^ vanishes simultaneously for all 7^ |m| < Mq for a finite number of values of y, we use 

y = — - ^ — , and check whether ([8| is well conditioned. If ^ is not well conditioned, we reduce y slightly. 

The algorithm ensures that we never reduce y by a factor of 9/10 or more. Now we check whether ^ vanishes 
simultaneously for all 7^ \m\ < Mq for the given y. If ^ does vanish, we continue with a finite number of random 

choices for the value of y € ( — — - — — - — , min Im w and check for each value of y whether ((8|) vanishes for all 



7^ |m| < M, 



10 2-kMq 



2.7. Turing's method. Turing's method is a method of verification that the list of eigenvalues of Maass cusp forms 
is consecutive, once we have a suitable bound for the error term S(t) in "average" Weyl's law 3\f(<) = M(t) + S(t). 
Roughly speaking, the method is the following. Assume that the error term S(t) in the "average" Weyl's law for 
the corresponding surface satisfies a bound of the type 

T 

E t (T) < (S(T)) := iy S(t)dt < E U {T), 


where Ei(T) —> and E U (T) — > 0, as T — > 00. Then, we have the following test of consecutiveness [51 13"3]. 

Step 1. Compute N num (T); the number of numerically found eigenvalues in the interval l/4<A<T 2 + l/4 and 
denote by 

r m (T) := K num (T) - M(T) 

the difference between the number of numerically found eigenvalues and the average Weyl's law. 

Step 2. Add a "fake" eigenvalue Af a k e near the end of the list of eigenvalues and compute (5 num (T)). If the value 
(S num (T)) exceeds E U (T), then the list of eigenvalues is consecutive in the interval 1/4 < A < Af a k c - 

3. Average Weyl's law for r (A r ) + 
In this section we prove Theorem [T] 

Let us recall that N — p\ ■ ■ ■ p r is a square-free positive integer, and define the function 

un{j,T) := Tlogpj - [ Jtt, 
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where, as previously stated, [x\ denotes the greatest integer less than or equal to x. Let Xm — ro(iV) + \H be the 
Riemann surface associated to the Fuchsian group Iq (N), and let Zx N denote the Selberg zeta function associated 
to X N . 

Let A <E (1, 3/2) and T > 1 be arbitrary real numbers, and let R(A) be the rectangle with vertices 1 — A — iT, 
A — iT, A + iT, 1 — A + iT. Without loss of generality, we assume that A and T are such that Zx N (s) ^ for 
s € dR(A). Formula (5.3) on p. 498 of [17] states the location of zeros and poles of the Selberg zeta function 
Zx N - In the notation of [17], one has that m = and W = id. Furthermore, <^jv(1/2) = —1, hence, application 
of Theorem 4.1 on p. 482 and formula (4.6) on p. 485 of [T7] yields that, in the notation of formula (5.3), one has 
A + B - K - C = 2. 9 - 2. Therefore, 

1 f Z' 

— / - 7 ^ L {s)ds = 2Ji N [Q < r n <T]+ 2Q N [0 < lm(p) <T]+2g-2 + n N + 2m 1/4 N , 

27Tl JdR(A) Z X N 

where Qat[0 < Im(p) < T] denotes the number of zeros p of the scattering determinant ipjq with Im(p) G (0, T]. 

Let dP(A) denote the polygonal path joining points 1/2 — iT, A — iT, A + iT and 1/2 + iT. Using the functional 

z' 

equation for the function 2)jv(s) := „ Xn (s), as in the proof of Theorem 2.28 on pp. 466-467 of 1171 . we can write 

Z ' X N 

(9) NN(T)+Q N [Q<lm(p)<T] = l-g- 1 ^-m 1/4 , N + R 1 (T) + - i - f ^(s)ds - / G N (s)ds, 

2 4m J ifN 4m 



dP(A) dP{A) 



where 



and 



i?i(T):=^. J V N (s)ds 

dP(A) 



(10) C„0O = Vol(X N )(s - 1/2) tan(.( S - 1/2)) - £ if * ~ " 1/2) 



t=i j 



-i mi sm(itj/rrii) C0S7r(s — 1/2) 



+ 21og2+^(l/2 + S ) + ^(3/2- S ). 



By Theorem 2.29 on p. 468 of [T7], we have the estimates 



T 



(11) ^ 1 ^ = (iogr) and J Rl ^) dt = asT ^°°- 

g q og 

To see that our function Ri(T) is equal to the function S(T) in Theorem 2.29 of [T7], we refer to Definition 2.27 on 
page 465 of [T7] . In addition, one can easily prove that Si (T) , in the notation of [T7] , coincides with the integral of 
Ri(T). To do so, one simply integrates the formula for i?i(T), interchanges the order of integration, evaluates the 
inside integral, and then integrates by parts. We choose to omit the details of these elementary calculations. 

Since A <= (1,3/2), the function Cjv(s) has no poles on the sides of the rectangle Ri/2(A) which has vertices at 
the points 1/2 — iT, A — iT, A + iT and 1/2 + iT. Furthermore, the only pole of Cat(s) inside Ri/2(A) is a simple 
pole at s = 1. Since 

lim — = 

a-tl COS7r(s — 1/2) 7T 

we conclude that 

Res^C^)- L rn^jJrnA 2^ + ^ I, 1 " m J ~ 1 " ^ 

Therefore, by the calculus of residues, having in mind that Cjv(l/2 + it) = 6^(1/2 — it), for real and non-negative 
t we get 

T 

(12) i j e N (s)ds = ~(l-2g) + ~ J e N (l/2 + it)dt. 

dP(A) 
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By substituting s = 1/2 + it into (10), we then have that 

T T T 

(13) fe N (l/2 + it)dt=-Vol(X N ) [ttwh(nt)dt-Y — ^ , l - . / C0Sh ^ (2j / m '" 1)t ^ 
J J TTii sm(n j/rrii) J cosh7rf 

+ 2 Re ^ ( 1 + it)dt j + 2T log 2 = 7 X (T) - 7 2 (T) + 7 3 (T) + 2T log 2, 

where, in obvious notation, 7i, 7 2 and 7 3 are defined to be the integrals in (13 1. We will now estimate each of these 
integrals. 

We write t tanh(7rf) = t— 2i/(l + cxp(27rf)) to get the expression 

/ oo \ 

h(T) = Vol(X N ) ^ - 2 J + 2 91 (T) 



where 



(14) 9i{T)=J- 



tdt 

_i_ plixi ' 



Quoting formula 3.411.3 from [IB] with v = 2 and /i = 27r, having in mind that ("(2) = 7r 2 /6 and T(2) = 1, we get 
(15) h(T) = -Vol(X N ) - i + 2 fll (T) 

Similarly, by quoting formula 3.511.4 from |16j with a = ir(2j /rrii — 1) and & = 7T, we arrive at the equation 

T 

cosh 7 r(2j/m l - l)t 1 . 



where 



Hence, 



cos!i7rf 2sin(7rj/mi) 



, . , . cosh 7 i-(2j/m l - l)i 
ff2( ^' T) := / dL 



i=i 1 



where we define 



02 cn •= v — v g ^ 

2 TOi j=i 8in ( 7r i/ m i) 

Finally quoting formula 8.344 from |16| . which is essentially Stirling's formula, with z — 1 + iT and n = 2 we 
get that 



(17) i7 3 (T)=Ref-i J Qogr(l + it))' dtj = Ln(Iogr(l + *T)) = -T + TlogT + J 



53 (T), 



where 

(18) ft(T) = 1 Im (log (l + ^)) + Re (t1o S (l + ^)) - + Im(^(T)) 
and 

(19) \MT)\ < — — < ^ 



12(1 + T 2 ) 3 / 2 cos 3 (i arg(l + iT)) " 12V2Tcos 3 (^/4) 180T' 
for T > 1. In the above computations, B 2 = 1/6 and B4 = —1/30 are Bernoulli numbers. 
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(20) 



Substituting (15), (16) and (17) into ([T3|) 3 and in turn using (12), we get the expression 

riii — 1 

] V_ 

2 24 



1 



n / x , 1/, „ s VoI(Xat) /T 2 1 
e JV (a)ds=-(l-2 ff ^ 1 

Z Z7T 



I , m ; — 1 

^ 4m,- ^ 



1 



ap(A) 



1 



TlogT-T 



tt\ log 2^ Vol(Xjv) 



T- 



'5i(7 1 ) + -ff 2 (T) + - ff3 (r). 

7T 27T 7T 



Using the evaluation pj of the scattering determinant, we immediately deduce that, inside the rectangle Ri/ 2 (A) 
the function ipn(s) has a simple pole at s = 1 and zeros at points p. Therefore, 



(21) 



1 



lfN -{s)ds = Q N [0 < Im(p) < T] - ' ' ' 



dP(A) 



2 4tt J ifiN 



(1/2 + it)dt. 



Combining (|2lJ with Q and |9]) yields, for T > 1, 



(22) N 7V [0<r„<T]=i? 1 (T) 



1 



^ (1/2 + it)dt + MM T2 - ^ + ^(1 - log 2) 

47T 7T 7T 



mi — I 



47T J (^JV 
— T 

Vol(Xjv) 1 njy 



^ 4m, ^ sin^Trj/m,) 



48tt 



- ™1/4,V 



/Vol^v) m 1 m 1 m 
Si CO - —gi{T) - - 53 (T) 

V 7T 27T 7T 



Taking logarithmic derivative of pi), we get 



(23) 



1 



47T J </?V 



^ (1/2 + it)tft = 



T 

1 /" D' f 



Am ] D 



(1/2 + it)dt. 



We now will compute the three integrals on the right-hand side of (23) separately. First, 

T 



(24) 



1 



dt 



1 1 



2tt J (1/4) + t 2 2 7T 



arctan(l/T). 



As for the second term on the right-hand side of ( 23 1 , we begin by writing 



-— Re [ —i [ (logf(l + 2i*))'dt I = -- Im(log£(l + 2*T) - log£(l)). 



From the definition of the function £, one has £(1) = £(0) = 1/2, so then 

(25) -- Re f -i [ (log£(l + 2it))'dt = -- Im(log£(l + 2iT) - log£(l)) 



-1 - + "(1 + logTr) - -g 4 (T) - -R 2 (T), 

TT TT 7T 7T 



where H 2 (T) = Im(logC(l + 2iT)) = Im(log(2iTC(l + 2iT))) - tt/2. From Stirling's formula, we have that 



(26) 



<7 4 (T) = Im log 1 



2iT 



Re Tlog 1 



J_\\ _ 2B 2 T 
2iT J J (1 + 4T 2 ) 



Im(K 4 (T)). 



The error term 3? 4 (T) satisfies the inequality 
(27) |X4(T)| < 



|S 4 | 



< 



12(1/4 + T 2 ) 3 / 2 cos 3 (iarg(l/2 + iT)) " 225T' 



for all T > 1, which we have deduced in a manner similar to (|19|) 
by 



(28) 



As for the third integral in (23), we begin by noting that the logarithmic derivative of the function Dn is given 

fa - l)logj>j ■/>: 



D N 



{ S ) = -\ogN-Y J 
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Furthermore, straightforward computations yield the formula 



(29) 
where 



1 

2^ 



Re 



(Pj - l )P 



1/2+it 



(Pj +P]){Pj +1) 



dt 



2?r log pjipj i + l) 



Tlog Pj 



1 + aj cos u 



a, = 2/(VP7+ l/VPj) = l/cosh((l/2)log(p,))- 



With these preliminary computations, the third term on the right-hand side of (23) can be evaluated using (28) 
and ( 29 ) , namely we have the formula 

(30) 



\QgN T 1 

2tt 2tt ^ pi + 1 

.7 = 1 



Tlogp 3 



f/ii 



1 + cos u 



We write 



(31) 



Tlogp 3 



1 + dj cos u 



£ . I 



fc=0 



_ , I T log Pj . 



+ dj COS It 



r/i/ 



T log Pj 

1 + ( — 1) L ir J ttj COS It 



and use [IB] , formulas 3.613.1 with n = 0, a = q ? and 2.553.3 with a = 1, 6 = (— 1)L s J Jaj (hence b 2 < a 2 ) to 
evaluate the two integrals in (31). Substituting (31) into (30), and employing the definition aN{j,T) := Tlogpj — 
\ Ty °^ V3 J 7r, we get the expression 



_L|£k (1/2+jt)(ii 

-T 



log TV 1 -A , . 1 <A 

-T + — a N {j, T) arctan 



7T 2?T 



3=1 



3=1 



JjJ-1 



/Pj + l 



(-1) 



tan 



aN{j,T) 



Now, by combining this last formula with (23), (24) and (25), we arrive at the expression 

- J ska, 



(32) 



4tt 7 </3jv 

-T 



(1/2 + it)dt 



1 TlogT T 



7T 7T 



(1 + log(7r/iV)) - -R 2 (T) + i- T) 

7T Z7T » — ' 



3=1 



— > arctan 

7T 



3=1 



1 

/p7+i 



(-i) L — ^ — J 



tan 



/ a N (j,T) 
V 2 



- (arctan(l/T) + gJT)) . 

IT 



Substituting ( 32 ) into ( 22 ) , we immediately see that 

>hr[0 < r n < T] - Mjv(T) = Sn(T), 

where 

M„(T) = ^&>r> - + I (2 + ,„ gW2W)) + £ -L "f 1 i 

v y 47T 7r 7r v v ' " ^ Am, ^ sin 2 f7n 



V6L(X N ) 



^4m, ^-J sin (nj/mi) 48tt 



l~ ^ + ^l^ aN ^ T > ~ arCtan 

3=1 3=1 



/p7 + i 



tan 



/ ajy(j, T) 
V 2 



TOi/ 4 ,JV 



+ Gat(T), 



with 

and 
(33) 



Gat(T) = (-2 Vol(Xjv) ff i(T) + fl2 (T) + 2 ff3 (T) + 2 ff4 (T) + 2 arctan(l/T)) 

Z7T 



Sn(T) = Ri{T) - -(Im(log(2iTC(l + 2iT))) - tt/2). 

7T 
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At this time, it remains to derive bounds for the error terms Gn(T) and Sn{T). From the definition (14 1 of the 
function g\ (T) we deduce that 



\ <h (T)\, lte-™dt= e -^-{T+^) 



For an arbitrary positive constant A > 0, the function f(x) = x 2 exp(^4 — Ax) is decreasing for x > 2/ A; hence, if 
A > 2, then f(x) < /(l) = 1 for all x > 1. Therefore, for A > 2, one gets exp(— Ax) < exp(— A)x~ 2 for all x > 1. 
Taking ^4 = 2ir > 2, we obtain the bound 



(34) 



e" 2 ^, 1 . 2tt + 1 1 



for all T > 1. Since uexp(l — u) < 1 for all « > 0, we get the inequalities 



\92(i,j,T)\ < - / exp((|(2j/m < ) - 1| - 1) u)du 



2 exp((|(2j/m i )-l|-l)7rr) 

7T (l-|(2i/mi)-l|) " (l-|(2j7m 4 )-l|) 2 7r 2 e 2" 



< 



Z mi — 1 

P77 - 



i 2e7r ^-^ sin(7rj/m.j) T 



For j € {1, . . . , mi — 1} one has 1 — |(2j/?7ij) — 1| > 2/m.i, hence 
i mj-l „ 

(35) | 52 (r)i < y — y 2 — 

S m i ^ e- sin^/m*) (1 - \(2j/ mi ) - l\f nT 
for all T > 1. 

In order to obtain bounds for g% and 34, we need to estimate Im (log (l + ^)) and Re (Tlog (l + ^j) for a = 1 
and a = 1/2. When T > 1 one has |a/iT| < 1, so then 

(— l) fc_1 / a 



k=l 



Therefore, 



Similarly, 



00 / *\t. 

to ( log ( 1 + if))l = S^TTT 



fc=l 



(-l) fc / . 1 

T 



< 



T' 



B.(rv«(i + ±)) 



(-l) fe - x /a\2fe 



fc=i 



2/c VT 



< 



2T' 



Now, from (|18j), (|19j), (|26j) and (|27j) we conclude that for T > 1 

1 1 



(36) 
and 
(37) 



\9s(T)\ < 
MT)\ < 



1 



1 



49 



2T 2T 12T 180T 45T 

1 1 1 2 _ 1291 

2T + 8T + 12T + 225T ~~ 1800T 



Finally, for T > 1 one has arctan(l/T) < 1/T, hence substituting (34 1, (35 1, (36 1 and (37) into the definition of 
Gat(T'), we arrive at 



\Gn(T)\ < 



1 / Vol(XAr)(27r+ 1) 



1 



5051 



2irT \ 27r 2 exp(27r) 



l nij-l 

E rrii \ - 
2^ ^ sinfTrj/mi) ' 900 I ' 



which is the inequality stated in ([I]). 

By (11 1, the proof of the theorem will be complete once we show that 



(Im(log(2^C(l + 2it))) - ir/2)dt = O 



T 



log 2 T 



as T-> 00. 
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In fact, we will prove the stronger bound 

T 



(38) J (Im(log(2itC(l + 2i*))) - ir/2)dt = O(logT) as T -> oo. 



o 

By the changes of variables s = 1 + 2it, we can write 

T l+2iT 

log(2itC(l + 2ii))dt = — [ log((s - l)C(s))ds. 
2i J 
l 

The function log((s — l)C(s)) is holomorphic in the closed rectangle with vertices 1, A, A+ 2iT and 1 + 2iT, so, by 
Cauchy's theorem, we have that 

T AT 

(39) { log(2<(l + 2ii))dt = j. j log((cr - l)C(a))da + j \og{{A - 1 + 2it)Q{A + 2it))dt 



2i 

It remains to estimate Im(Ji(T)) and Im(J2(T)). Trivially, one has 



~ J log((CT + 2iT - 1)C(ct + 2iT))da = 0(1) + J X (T) + J 2 (T) as T -> oo. 



1 log (jK* (l + H + Im I y log ((A + 2it) 

It is elementary to show that Im (log (2it (l + ^r))) = n/2 + 0(t~ 2 ) for t ^> 1. The Dirichlet series representation 
of logC(^4 + 2zt), is absolutely and uniformly convergent in the range under consideration since A > 1. Therefore, 
we get the bounds 

(40) Im(J 1 (T)) = ^T + Q(l)+^ W 2 Im ( — ^ )=£r+0(l) as T — > oo, 

2 *— * n A log n V 2i J 2 



where A(n) is the von Mangoldt function. Combining (40) with (39), we have that 
T 

/"(Im(log(2itC(l + 2ii)))-7r/2)di = Im(J 2 (T)) + 0(l) as T oo. 
o 

In order to prove p8| ), we need to show that Im(J 2 (T)) = O(logT) as T — > oo. The proof of this bound is 
straightforward. Simply combine the elementary bound log(tr + 2iT — 1) = O(logT) together with the estimate 
log£(er + 2iT) = O(logT) which holds uniformly for a € which we quote from Theorem 3.5 in [32] . 

With all this, the proof of Theorem [l] is complete. 

Remark 14. In order to prove Corollary [2j one follows the analysis above up to equation (33). At that point, one 
uses the first part of Theorem 2.29 on page 468 of [T7] to bound the first term and Theorem 3.5 in to bound 
the second term. 

We now state three special cases of Theorem [l] first when T = PSL(2,Z), next when T = To(5) + , and finally 
when T = r (6) + . 

Corollary 15 (Average Weyl's law for PSL(2,Z)). 

N 1 [0<r„<T]-M 1 (T) = 5 1 (T), 

where 



with 



M 1 (T) = j^T 2 + T + log(7r/2)) _ m + n 

iz IT IT 144 



1 / 2tt + 1 + 6(1 + 2^3) exp(27r-l) 5051 \ 11 
1 11 n ~ 2tt \ 67rexp(27r) 900 j T T 
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and 



J Si(t)dt = O 



T 



\oe 2 T 



, as T — > oo. 



Proof. We apply Theorem [I] with N = 1. In this case, Di(s) = 1 and the signature of the group is (0;2,3;1). 
Furthermore, Ai > 1/4, by Theorem 11.4 from [21], hence rii = 1 and miu t jy = 0. □ 

Remark 16. We have been informed that in [8j, the authors prove an average Weyl's law for SL(2,Z) together 
with effective bounds for the integral of Si , using a trace formula approach. 



Corollary 17 (Average Weyl's law for r (5)+). Let a 5 (T) = Tlog5 



Tlog5 



7T. Then, 



wht 



N 5 [0<r n <T]-M 5 (T) = 5 5 (T), 
T 2 2TlogT T,_ , /7T\, 43 a 5 (T) 



M5 ^ = T ^ + ^ 2 + 1 ° g (lo) ) -48 + 



2?r 



1 



arctan 



'5 + 1 



tan 



« 5 (T) 




G 5 (T), 



and 



|G 5 (T)| < 



1 /27r + l + 6exp(27r-l) 5051 A 1 1 



2tt 



27r exp(27r) 



900 7 T T 



S 5 (*)di = O 



T 



log 2 T 



, as T->oo. 



Proof. We apply Theorem [T] with iV = 5. In this case, the signature of the group is (0;2,2,2;1). Also, the only 
eigenvalue < 1/4 is Ao = 0, by Lemma 10 and Corollary 11.5 from [21]. □ 





Tlog3 


7T 




7T 





7T. The 



where 
M 6 (T) 



N 6 [0<r n <T]-M 6 (T) = 5 6 (T), 



1 

arctan 

7T 



arctan 




(-i) L 



(~i) L 



tan - ( Tlog2- 



tan- ( Tlog3 



Tlog2 



Tlog3 



with 



and 



\G 6 (T)\ < 



1 /27r + l + 6cxp(27r~l) 5051 A 11 



2tt 



27r exp(27r) 



900 7 T T 



S 6 (t)dt = O 



T 



log 2 T 



, as T->oo. 



Proof. The proof is a straightforward corollary of Theorem [T] and basic properties of Fq(6) + . 



□ 
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n 


A„ for r (5) + 


A„ for r (6) + 


1 


17.32676 


20.93844 


2 


24.23291 


26.24717 


3 


36.89998 


37.71537 


4 


40.58784 


40.01593 


5 


46.81219 


52.39092 


3555 


15623.315 


15649.988 


3556 


15623.860 


15654.937 


3557 


15625.094 


15665.201 


12470 




52875.046 


12471 




52876.076 


12472 




52879.257 


12473 




52894.324 


12474 




52899.011 



Table 1. Eigenvalues of the Maass cusp forms on ro(5) + in the interval < A < 125 2 + 1/4 and 
on r (6)+ in the interval < A < 230 2 + 1/4. 



4. Numerical computations 

In this section we present numerical results on computations and statistical distribution of large sets of consecutive 
eigenvalues of Maass cusp forms on X$ and X§ . 

4.1. Computation of consecutive list of eigenvalues of Maass forms on X§ and Xq. A systematic search [31] 
for Maass cusp forms on r (5) + in the interval < A < 125 2 + 1/4 and on r (6)+ in the interval < A < 230 2 + 1/4 
results in 3557 and 12474 Maass forms, respectively. A few eigenvalues are listed in Table [T] At some point, the 
entire list of eigenvalues will be made publicly available. Prior to that time, the list will be made available to anyone 
upon request. 

We note that the lowest point of the fundamental domain of the surafce Xg has a larger imaginary part than 
that for A5. The height y of the lowest point has an influence on how many terms are to be considered in the 
Fourier expansion This is the reason, why the computations were much faster on X e than on X 5 . 

The algorithm for computing eigenvalues is described in detail in |3T] . The main ingredients are the following. 
First, using a set of trial values Ai, . . . , A„, we linearize Hejhal's system of equations ^ in the eigenvalue A around 
each trial value A. For each A, we obtain a matrix eigenvalue equation which is then solved numerically. In this 
step, the eigenvalues A of Maass cusp forms are related to the matrix eigenvalues via perturbation theory. As a 
result, one obtains a preliminary list of potential eigenvalues of Maass cusp forms. For each potential eigenvalue, 
we solve Q for 7^ \m\ < Mq and check whether the corresponding non-trivial solution is indeed a Maass cusp 
form. 



The check is described in section 2.6 This results in a verified list of Maass cusp forms. Finally, we need to 
check and verify that the list of Maass cusp forms is consecutive. As stated, this check is performed using "average" 
Weyl's law and Turing's method. If it turns out that eigenvalues are missing, we search for them, using additional 
trial values A, until our list of Maass cusp forms becomes consecutive, as indicated by Turing's method. 

We do not have rigorous Turing bounds, yet. Therefore, we use Turing's method heuristically. In light of the data 
obtained, and presented in various figures in this section, let us again discuss Turing's method, this time keeping 
the figures in mind. 

Let 3\f™ m (T) count the number of numerically found eigenvalues in the interval 1/4 < A < T 2 + 1/4. The 
difference between the number of numerically found eigenvalues and the average Weyl's law 

S» um (T) := W™ m (T) - M N {T) 



is a fluctuating function. Its mean comes close to a non-positive integer whose absolute value counts the number of 
solutions which have been overlooked. 
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Figure 5. The mean (S£ um (T)). 
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FIGURE 6. Mean (S£ um (T)), with the eigenvalue A 936 7 = 200.0359 2 + 1/4 removed. 
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Figure 7. Mean (Sg um (T)), with the fake "eigenvalue" A = 200 2 + 1/4 inserted. 



Figures [2] and [3] show the fluctuations S^ um (T). In figures g and [5j the mean 



(S^(T)) := - / S^(t)dt 



tends to zero for large T which indicates that all solutions have been found numerically. 

If a solution would have been overlooked, the graph would deviate from zero quite significantly. A demonstration 
is given in figure [6j where we have intentionally removed the eigenvalue A9367 = 200.0359 2 + 1/4, whereas in figure 
[7J we have intentionally inserted a fake "eigenvalue" at A = 200 2 + 1/4. 

T 

If we would have an explicit and efficient upper bound on J Sjsr(t)dt, we could apply Turing's method to prove, 



not just verify, that the numerically found lists of eigenvalues are consecutive. The proof would be to add a fake 
"eigenvalue" near the end of each list of eigenvalues and show that with this extra "eigenvalue" (S^ m (T)) would 
exceed the upper bound, as explained in section 2.7 see also |331 15115]. In our notation, what is needed is to explicitly 



evaluate the implied constant in the average of Sn(T). The algorithm in [T3] and [55] does, in fact provide such a 
bound, but the explicit value is somewhat large, hence impractical. 



Remark 19. For computing S'™ m (T) we need to evaluate Mjv(T) which includes the term Gjv(T). Actually, we 
do not know the exact value of Gn{T). According to the bound ([!]), we can safely neglect Gn{T) in the evaluation 
of S^ m (T) for T large. 

By Theorem]!] Mat(T) includes terms which depend on ajs[(j,T) and on arctan(-). For evaluating the average 
(Sif m (T)), we need to integrate over these terms. The sum of the ajy(j, T) and the arctan(-) dependent terms is 
periodic. We perform the integration by expanding the periodic contribution into a Fourier series, integrate the 
individual Fourier terms, and then sum up numerically. 

4.2. Nearest neighbour spacing statistics. Concerning the statistical properties of the eigenvalues, we must 
emphasize that the conjectured properties depend on the choice of the surface X. Depending on whether the 
corresponding classical system of a point particle that moves freely on the surface is integrable or not, there are 
some generally accepted conjectures about the nearest neighbour spacing distributions of the eigenvalues in the 
limit A — > 00. 

Whenever we examine the distribution of the eigenvalues we consider the values on the scale of the mean level 
spacings. 
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Figure 8. Nearest neighbour spacing distributions P(s) for the moonshine groups To(5) + (left), 
and To(6) + (right), which come close the Poisson distribution Pp isson(s) = e~ s . 

Conjecture 20 ( 2 ). If the corresponding classical system is integrable, the eigenvalues behave like independent 
random variables and the distribution of the nearest neighbour spacings is in the limit X — > oo close to a Poisson 
distribution, i.e. there is no level repulsion. 

Conjecture 21 ([4 ]). If the corresponding classical system is chaotic, the eigenvalues are distributed like the 
eigenvalues of hermitian random matrices. The corresponding ensembles depend only on the symmetries of the 
system: 

• For chaotic systems without time-reversal invariance the distribution of the eigenvalues approaches in the 
limit X — > oo the distribution of the Gaussian Unitary Ensemble (GUE) which is characterised by a quadratic 
level repulsion. 

• For chaotic systems with time-reversal invariance and integer spin the distribution of the eigenvalues ap- 
proaches in the limit X — > oo the distribution of the Gaussian Orthogonal Ensemble ( GOE) which is char- 
acterised by a linear level repulsion. 

• For chaotic systems with time-reversal invariance and half-integer spin the distribution of the eigenvalues 
approaches in the limit X — > oo the distribution of the Gaussian Symplectic Ensemble ( GSE) which is 
characterised by a quartic level repulsion. 

These conjectures are very well confirmed by numerical calculations, but several exceptions are known. 

Exception 22. The harmonic oscillator is classically integrable, but its spectrum is equidistant. 

Exception 23. The geodesic motion on surfaces with constant negative curvature provides a prime example for 
classical chaos. In some cases, however, the nearest neighbour distribution of the eigenvalues of the Laplacian on 
these surfaces appears to be Poissonian. 

With our lists of consecutive eigenvalues, we can examine the nearest neighbour spacings. We unfold the spectrum 

u n = M N (r n ) with A„ = r 2 n + 1/4, 

in order to obtain rescaled eigenvalues u n with a unit mean density. Then 

defines the sequence of nearest neighbour level spacings which has a mean value of 1 as n — > oo. For the moonshine 
groups ro(5) + and To(6) + we find that the spacing distributions come close to that of a Poisson random process, 

-^Poisson (^) 6 , 

see figure |8j as opposed to that of a Gaussian orthogonal ensemble of random matrix theory, 

Pqoe(s) ^ | se ~ fs2 - 

The spacing distributions are in accordance with Conjecture [3] 

One might wonder whether eigenvalue spacings are correlated. For this we investigated joint eigenvalue spacing 
distributions, 

P(s,s')dsds' = P((s n ,s n+1 ) e [s,s + ds) x [s 1 , s' + ds')). 
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FIGURE 9. Joint nearest neighbour spacing distribution P(s,s') for ro(5) + which comes close to 
a product of Poisson distributions Ppoisson(s)Pp isson(s')- 




Figure 10. Joint nearest neighbour spacing distribution P(s,s') for To(6) + which comes close to 
the product of Poisson distributions -Pp isson(.s)-Ppoisson(s')- 

For ro(5) + and ro(6) + , we find that the joint eigenvalue spacing distributions factor into a product of Poisson 
distributions, 

P{s,s) = Ppoi SS on(s)fpoi SS on(s'), 

see figures [9] and [TU] Hcuristically, the spacings between rescaled eigenvalues are uncorrelated which implies that 
the eigenvalues are uncorrelated as well. 
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5. Concluding remarks 

5.1. Topological equivalence versus Weyl's law. The groups To(5) + and ro(6) + , which were the main focus 
of investigation in our paper, are topologically equivalent, have different Weyl's law, yet the two sets of eigenvalues 
seem to have the same spacing distributions. 

From the tables presented in [12] . one can find other examples of such groups. In the case when the genus g is 
zero, we have the following examples of topologically equivalent groups with different "classical" Weyl's laws and 
"average" Weyl's law. 

(1) For N G {11, 14, 15}, the signature of the surface r (iV)+\H is (0; 2, 2, 2, 2; 1), 

(2) For N e {17, 22, 30}, the signature of the surface r (A r )+\H is (0; 2, 2, 2, 2, 2; 1), 

(3) For N e {23, 33, 42}, the signature of the surface r (7V)+\H is (0; 2, 2, 2, 2, 2, 2; 1), 

(4) For N £ {29, 38}, the signature of the surface r (7V)+\H is (0; 2, 2, 2, 2, 2, 2, 2; 1), 

(5) For N £ {46, 51, 55, 66, 70}, the signature of the surface r (iV)+\H is (0; 2, 2, 2, 2, 2, 2, 2, 2; 1). 

There are also examples of the groups with genus g — 1, such as N <E {83, 123, 143, 182, 195} each of which has 
signature (1; 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2; 1) and whose Weyl's law asymptotics differ in the T term. 

This empirical investigation yields to an interesting question: For a given positive integer k, is it possible to find 
k topologically equivalent surfaces arising from moonshine groups having different Weyl's laws? 



5.2. Weyl asymptotics versus nearest neighbour statistics. Generally speaking, discrete eigenvalues of the 
Laplacian, or, equivalently, positive imaginary parts of zeros of the corresponding Selberg zeta function on the critical 
line, are increasing sequences of numbers, and the associated Weyl's law is an approximate counting function of 
such sequences. The results in section [4. 2| are re lated to numerical computation of the nearest neighbour statistics 
of eigenvalues of Maass cusp forms on ro(-/V)+\EI, for TV — 5 and N ~ 6. We have seen, empirically, that the 
nearest neighbour statistics for the eigenvalues of Maass cusp forms seem to be each equal even though the Weyl's 
laws are different. One may argue that the reason for this is that the Weyl's law differs in the T term, while the 
first two lead terms are the same in the two cases we considered. 

Therefore, a natural question which arises is to what extent does the nearest neighbour statistics of an increasing 
sequences of numbers depend on its average counting function. The answer to this question is presented in the 
following example. 

Example 24. Let {x„}„ e N be an increasing sequence of numbers having a mean density of 1, by which we mean 

lim i#{x„ < T} = 1. 

T— >oo 1 

Let m(t) be an increasing function, defined for t > 0, such that m(0) > 1/2. (In the Weyl's law case, m(t) = 
dot 2 + aitlogt + a^t + 03 + . .., for some positive number do-) Let us define a sequence of numbers A„ by letting 
A„ := m~ 1 {x n — I), where mT 1 denotes the inverse function of m. Let be the counting function 

>T(t) := #{A„ < t}, 

and let the Weyl asymptotics M(f) be a smooth approximation to N(t) such that 

T 

lim i / CN(t) - M(t))dt = 0. 

T—foo 1 J 


The unfolded spectrum {u n } is defined by u n := m(A„). Trivially, u n = m(A n ) = m(m~ 1 (x n — |)) = x n — ^, for 
all n € N hence u n+ i — u n — x n+ i — x n , so the nearest neighbour statistics of the unfolded spectrum {u n } equals 
the nearest neighbour statistics of the initial sequence {x n }. 

We are free to distribute the sequence of increasing numbers {x n } such that the nearest neighbour statistics 
of {x n } coincides with our favorite distribution of non- negative numbers. We are also free to choose the smooth 
increasing function m(t), and hence the Weyl asymptotics arbitrarily. Since {x n } and m(t) can be chosen indepen- 
dently of each other, we conclude that the nearest neighbour statistics of the unfolded spectrum {u n } is completely 
independent of the Weyl asymptotics. 

Therefore, all the analytic results on the Weyl asymptotics are completely independent of the numerical results 
on the nearest neighbour statistics. Neither carries any information of the other, regardless of how many expansion 
terms we include in the Weyl asymptotics. Analytics and numerics complement each other. 
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